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We report first-principles density-functional pseudopotential calculations on the atomic structures, 
electronic properties, and band offsets of BaO/BaTi03 and SrO/SrTiOs nanosized heterojunctions 
grown on top of a silicon substrate. The density of states at the junction does not reveal any elec- 
tronic induced interface states. A dominant perovskite character is found at the interface layer. The 
tunability of the band offset with the strain conditions imposed by the substrate is studied. Using 
previously reported theoretical data available for Si/SrO, Si/BaO and BaTiOs/SrRuOa interfaces 
we extrapolate a value for the band alignments along the whole gate stacks of technological interest: 
Si/SrO/SrTi0 3 and Si/BaO/BaTi0 3 /SrRu0 3 heterostructures. 
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I. INTRODUCTION 

The search for alternative gate dielectric materials to 
replace silica (Si02) in microelectronic devices is one of 
the grand challenges that the materials science commu- 
nity and the Si-based semiconductor industry are facing 
at the current time H. The rapid scaling of the physical 
gate lengths of Metal-Oxide-Semiconductor Field-Effect- 
Transistors (MOSFET) requires a concomitant rapid re- 
duction of the gate dielectric thickness in order to pre- 
serve a high gate oxide capacitance. This can no more 
be accomplished by lowering the size of the SiC>2 layer 
because, together with problems in the thickness control, 
the leakage current would become inacceptably high. In- 
deed, the leakage current from the channel to the gate is 
due to the direct tunneling of carriers and increases ex- 
ponentially with the decrease of both the thickness of the 
gate dielectric and the height of the electrostatic barrier 
for the electrons through the gate stack. The current 
roadmap projection (assessed by the the International 
Technology Roadmap for Semiconductors, ITRS □) im- 
poses the choice of an alternative gate dielectric with a 
good capacitance for a thick-enough layer and its full im- 
plementation into the production line by 2005. 

The properties that the new dielectric should meet are 
well established and have been reviewed recently by Wilk 
and Wallaces. They can be divided into fundamental ma- 
terial properties, and device processing and performance 
properties. Amongst the material properties, we can enu- 
merate (i) a higher dielectric constant than amorphous 
silica (Ksi0 2 = 3.9) in order to increase the capacitance 
without decreasing the thickness, (m) large band gaps 
and band offsets with Si to prevent tunneling currents, 
(Hi) a good thermodynamic stability in contact with the 
Si substrate, (iv) a good quality of the interface with the 
Si channel, which means a small number of electrical de- 
fects and a low midgap interface state density, and (v) 
film morphology avoiding the formation of polycristalline 



films and grain boundaries. Amongst the device proper- 
ties, we can cite (vi) a good compatibility with metallic 
gate electrodes, (vii) a compatiblity with the deposition 
mechanism during the fabrication process, (viii) reliabil- 
ity. 

Many materials satisfy some subset of the previous cri- 
teria, but the identification of a dielectric that addresses 
simultaneously all of the requirements is a real challenge. 
Investigations on oxides like AI2O3, ZrC>2, HTO2, Ta^Og, 
Y2O3, Gd2C>3, and TiO*, have thrown encouraging re- 
sults in the last few yearsQ. Amongst the most promising 
candidates, ABO3 perovskite oxides (where A stands for 
Ba or Sr and B stands for Ti) appear in good position. 

The ABO3 compounds have a dielectric constant above 
300, one order of magnitude higher than the other can- 
didates. Although they are thermodynamically unsta- 
ble in direct contact with Si (they react to form tita- 
nium silicide and alkaline-earth silicateoO), they can be 
grown in perfect registry with the Si substrate by means 
of Molecular Beam Epitaxy (MBE) when including a 
silicon-compatible buffer layer. On one hand, this layer 
must be sufficiently thick to ensure the physical separa- 
tion between the substrate and the perovskite. On the 
other hand, it must remain thin enough to keep the bene- 
fit of the high dielectric constant of the ABO3 compound 
(the capacitance of the lower-K buffer layer being in series 
with that of the perovskite). | — . | — . | — . 

In the McKee- Walker processBLjfl, the buffer consists 
in few atomic layers of AO alkaline-earth oxide that can 
eventually be alloyed during the growth in order to ac- 
comodate the lattice mismatch with Si. The growth of 
AO on Si includes the presence, at the interface, of a 
(sub) monolayer of ASi2 silicide so that the final structure 
corresponds to the sequence Si/ASi2/AO/AB03. The 
epitaxy is such that AB0 3 (001) || AO (001) || Si (001), 
and ABO 3 (110) || AO (100) || Si (100), i. e. the AB0 3 
atomic planes are rotated 45° around the (001) AO di- 
rection The epitaxial crystalline growth at the ox- 
ide/semiconductor interface avoids the formation of de- 
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fects and ensures the continuity of the dielectric displace- 
ment El. MBE techniques allow the control of the grow- 
ing sequence at the submonolayer level preventing grain- 
boundaries and providing a good quality interface and 
extremly smooth surface morphology. 

First attempts to make MOSFETs including per- 
ovskite oxides have been reported recently. Using a 110 
A-thick SrTiC>3 layer as the gate dielectric, Eisenbeiser 
et al Ej have fabricated a transistor that behaves com- 
parably to a 8 A-thick Si0 2 /Si MOSFET. The improve- 
ment in transistor performance was very satisfactory, and 
the leakage currents was two order of magnitudes smaller 
than in a similar SiC>2-based device. 

As it was pointed out before, the barrier height of the 
dielectric with respect to the Si substrate should be large 
enough to minimize carrier injection into the conduction 
band states. A large value of the Conduction Band Off- 
set, CBO, between Si and the gate dielectric is required, 
and typically materials with CBO smaller than 1.0 eV axe 
rejected for further applications. Robertson andnChen Ell, 
aligning the Charge Neutrality Levels (CNL)B of both 
semiconductors, have estimated the CBO for a Si/SrTi03 
interface to -0.14 eV (SrTi03 below, that is no barrier at 
all for the elections) in very good agreement with exper- 
imental results^. This prevents, in principle, the use of 
the titanate as the gate dielectric in electronic devices. 
However, the presence of the buffer alkaline-earth oxide 
in the heterostructure was missing in their approach. In 
this paper, we will show that in addition to providing 
a physical separation between Si and the perovskite, the 
presence of the alkaline-earth oxide also allows to monitor 
efficiently the band offset. 

We report a study of the properties of BaO/BaTi03 
(from now on, we will refer to this heterostructure as the 
Ba- interface) , and SrO/SrTi03 (Sr- interface) structures 
from first-principles. The method on which the simula- 
tions are based is described in Section ||. In Section 



[II, we discuss the details of the atomic structure at the 



interfaces. The electronic structure is presented in Sec- 



tions IV, where we analyze the density of states at the 
junctions. In Section ^ we stud y t he band-offset at the 
interface. Finally, in Section [Vl|, an estimate of the 
band alignment of the whole whole Si/SrO/SrTiOs/Pt 
and Si/BaO/BaTi03/SrRu03 structures will be given. 



core and valence states, the 3s and 3p electrons of Ti, 4s 
and Ap electrons of Sr, and 5s and bp electrons of Ba were 
explicitly included in the calculation. Ti, Sr and Ba pseu- 
dopotentials were generated scalar-relativistically. The 
reference configuration and cutoff radii for all the atoms 
we used are shown in Table-||. 

The one-electron Kohn-Sham eigenstates were ex- 
panded in a basis of strictly-localized E3 Numerical 
Atomic Orbitals IHJ. Basis functions were obtained by 
finding the eigenfunctions of the isolated atoms confined 
within the new soft-confinement spherical potential pro- 
posed in Ref. ^5|. We used single-^ basis set for the 
semicore states of Ti, Sr and Ba, and double-^ plus po- 
larization for the valence states of all the atoms. For Sr 
(respectively Ba) an extra shell of 4d (respectively 5d) 
orbitals was added. All the parameters that define the 
shape and the range of the basis functions for Ba, Ti and 
O were obtained by a variational optimization in cubic 
bulk BaTi03, following the procedure described in Ref. 
p5f For Sr, another optimization was performed in bulk 
SrTiOs, frozen in the atomic orbitaJs of Ti and O to these 
previously optimized in BaTi03 Ea. 

The electronic density, Hartree and exchange- 
correlation potentials, as well as the corresponding ma- 
trix elements between the basis-arbitals, were calculated 
in an uniform real-space grid An equivalent plane 
wave cutoff of 200 Ry was used to represent the charge 
density. Once self-consistency was reached, the grid was 
refined (reducing the distance between grid points by 
half) to compute the total energy, atomic forces and 
stress tensor. 

The integrals in reciprocal space were well converged, 
using in all the cases a sampling in k o£ihc same qual- 
ity as the (6x6x6) Monkhorst-Pack lA mesh in bulk 
BaTi03. The equivalent cutoff- length cB for this sam- 
pling, 13 A, was the one employed in all simulations. 
This represents a large number of fc-points thought that 
all the materials involved in the heterojunctions are in- 
sulators. However it has been proved that this fineness 
is mandatory while dealing with perovskites Ej. 

Test of the performance of the Siesta method on 
perovskites were done in bulk BaTi03t2l. Lattice con- 
stants, ferroelectric distorsions, Born effective charges, 
and phonon dispersion curv es are in very good agreement 
with plane-waves SBBe2 and full potential LAPW cal- 
culations a. 



II. TECHNICALITIES 



Our calculations have been performed within Den- 
sity Functional Theory (DFT) EJ and the Local Den- 
sity Approximation (LDA) EJ. We used a Numerical 
Atomic Orbital (jMAXIX method, as it is implemented in 
the Siesta code llaO'ta. The exchange-correlation func- 
nal was approximated using the Perdew,and Zunger 
parametrization of Ceperley- Alder dataE3. 
Core electrons were replaced by ab-initio, norm- 
conserving fully-separablecil Troullier-MartinE3 pseu- 
dopotentials. Due to the large overlap between the semi- 



III. ATOMIC STRUCTURE AT THE 
INTERFACE 

In Table |IJ we report the experimental and calculated 
lattice parameters of the different materials involved in 
our heterostructures, together with the lattice mismatch 
with respect to the Si substrate. The misfit is defined as 
/ = 100 x (a — asi) /asi, where a and asi are, respec- 
tively, the lattice constant of the epilayer and Si. The 
value of / is positive when the epilayer is compressed 
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TABLE I: Reference configuration and cutoff radii (in bohr) of the pseudopotentials used in our study. Because of the inclusion 
of the semi-core states in valence, and within the Troullier-Martin scheme, Ba, Ti, and Sr pseudopotentials must be generated 
for ionic configurations (ionic charge +2). However, these are more suitable than the neutral ones, given the oxidation numbers 
of these atoms in the alkaline-earth oxides and perovskites. 





Ba 


Sr 


Ti 


O 


Reference 


5s 2 ,5p 6 ,5-i°,4/° 


4 S 2 ,4p 6 ,4rf°,4/ 


3s 2 ,3p 6 ,3d 2 ,4/° 


2 S 2 ,2p 4 ,3d°,4/° 


Core radius (a.u.) s 


1.75 


1.50 


1.30 


1.15 


P 


2.00 


1.50 


1.30 


1.15 


d 


2.50 


2.00 


1.30 


1.15 


f 


2.50 


2.00 


2.00 


1.50 



TABLE II: Experimental and theoretical lattice constants (a, 
in A) for the different compounds involved in our heterostruc- 
tures. The lattice mismatch, /, between a given epilayer and 
the Si substrate (in % with respect to the substrate lattice 
constant) is also reported. cIa-a (A = Ba or Sr), stands for 
A-A nearest neighbour distance in AO oxides. Perovskite val- 
ues refer to the cubic structure. 



System 




Experimental 
a (A) / (%) 




LDA-DFT 

a (A) 


/ (%) 


Si 




5.43 a 




5.389 b 




BaO 




5.52 a 1.66 




5.433 


0.82 




(d Ba 


-Ba = 3.90) 


(d B a 


-Ba = 3.842) 




BaTiOs 




4.00 c 4.18 




3.948 


3.60 


Si 




5.43 a 




5.389 6 




SrO 




5.16 11 -4.97 




5.075 


-5.83 




(dsr 


-s r = 3.65) 


{dsr 


St = 3.588) 




SrTiOa 




3.91 d 1.83 




3.874 


0.90 



a N. W. Ashcroft, and N. D. Mermin, Ref. |i| 
'J. M. Soler et ai, Ref. EL 
C G. H. Kwei et al, Ref^M 
d T. Mitsui et ai, Ref. |37} 



and negative when it is expanded. In Table [II], we ob- 
serve that the LDA produces a systematic underestimate 
of the lattice constant (about 1%). Nevertheless, the cor- 
rect sequence of lattice mismatch is obtained so that the 
calculations will reproduce the experimental strain con- 
ditions when working at the theoretical lattice constants 
of the substrate. 

Interfaces were simulated using a supercell approxima- 
tion. The basic unit cell, periodically repeated in space 
corresponds to the generic (AO) n /(AO-B02) m formula, 
where n and m are respectively the number of ACLpxide 
atomic planes and the number of ABO3 unit cells 0. For 
even n and odd m (the only cases studied in this work), 
this structure possesses two mirror symmetry planes lo- 
cated on the central AO and BO2 layers. 

We considered pseudomorphic heterojunctions, so that 
the lattice constant parallel to the plane of the interface, 



TABLE III: Theoretical values of the elastic constants en 
and C12 in Mbar. 





en 


C12 


2(ci 2 /cii) 


BaO 


2.10 


0.57 


0.54 


BaTiOs 


3.71 


1.26 


0.68 


SrO 


2.36 


0.57 


0.48 


SrTiQ 3 


3.93 


1.17 


0.59 



an, is assumed to remain the same on both sides of the 
structure. The choice of an allows to treat implicitly the 
mechanical effect of the substrate, which is not included 
explicitely in the calculations. 

To establish the notation, we will call the plane parallel 
to the interface the (x, y) plane, whereas the perpendic- 
ular direction will be referred to as the z axis. 

Under the strain conditions imposed by the Si sub- 
strate, the epitaxial layers will minimize the elastic en- 
ergy by elongation or compression of the lattice constant 
along z, a±. To determine its value, strain relaxations 
of the bulk unit cells of AO and ABO3 were performed 
under the constraint of fixed an. Since the lattice misfit 
between the substrate and the epilayers is small enough 
to remain in the linear regime, the different values of a± 
with respect to the in-plane lattice constant can be pre- 
dicted from the Macroscopic Theory of Elasticity (MTE), 
and therefore an estimation of the atomic structure of the 
interface can be done. Following the description of Ref. 
|38| , and for an interface orientation along (001): 

«i,X = [l - A:£i,||] a% 



A = 2^ (l) 

where ai, c\ x and c\ 2 stand for, respectively, the equilib- 
rium lattice parameter and the elastic constants of ma- 
terial i. Theoretical values of the elastic constants are 
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TABLE IV: Lattice constant perpendicular to the plane of 
the interface, a± , at different values of the in-plane lattice con- 
stant, ay. Results from both, first-principles structural mini- 
mizations (FP) and macroscopic theory of elasticity {MTE) 
are reported. Units in A. 



System 


a || 


MTE 

a ± 


a ± 




O.ooy 


K A ^7 
0.40 / 


0.40 / 


i 


0.4OU 


0.4O0 


K A 

O.4oo 




5.665 


5.307 


5.322 




3.811 


4.041 


4.054 


BaTiOs 


3.839 


4.022 


4.025 




4.006 


3.909 


3.911 




5.389 


4.924 


4.939 


SrO 


5.430 


4.905 


4.923 




5.522 


4.861 


4.893 




3.811 


3.912 


3.915 


SrTi0 3 


3.839 


3.895 


3.893 




3.904 


3.856 


3.857 



constant was constrained to the theoretical one of Si. 
The most important features are: (i) a compression of 
the interplanar distance at the interface layer ; (ii) the 
appearance of an ionic interface dipole, due to the oppo- 
site motion of the anion and the cation at the AO-layer 
closest to the interface - the anion moves towards the 
AO region whereas the cations displaces inside the ABO3 
part -; (iii) a monotonic decay of the absolute value of 
the rumpling parameter as a function of the distance to 
the AO-interface layer, where the major relaxations are 
localized, and (iv) the oscillatory behaviour of the sign 
of r]i and (3i inside the perovskite from layer—taJayer, as 
it happens also in ABO3 free-standing slabsErEJ. 

The main difference between the Ba and the Sr- 
heterostuctures is the magnitude of the relaxations at the 
interface, larger in the last case. All these conclusions are 
independent of the in-plane lattice constant imposed in 
the calculation, and show very good-agreement with the 
results obtained using the Abinit e3 plane-wave pseu- 
dopotential code. 



reported in Table [II. Bulk structures from the macro- 



scopic theory are in excellent agreement with the first- 
principles results, as can be drawn from the results in 
Table IV (relative errors within 1 % for all the cases). 

The resulting bulk tetragonal unit cells were used as 
the building blocks of our supercell. However, as inter- 
planar distances in the region close,to the interface are 
not predicted properly from MTE cB, a full relaxation of 
the geometry using first principles methods was needed. 

For each interface, a reference ionic configuration was 
defined by piling up truncated bulk strained materials. 
Atomic coordinates were then relaxed until the maximum 
component of the force on any atom was smaller than 10 
me V/A. The maximum component of the stress tensor 
along z was smaller than 5 x 10~ 3 eV/A 3 for the Ba- 
interface and than 7 x 10~ 3 eV/A 3 for the Sr-interface. 
It has been confirmed that additional relaxation of a± 
(neglected in this work) does not produce any significant 
change. 

In order to characterize the atomic displacements in- 
duced by the relaxation, we define 5 Z (Mi) (respectively 
$z (Oi)) as the displacement of the cation (respectively 
oxygen) along z at layer i, with respect to the ini- 
tial reference configuration. We introduce the displace- 
ment of the mean position of each atomic plane as 
Pi = [$z (Mi) + S z (Oi)] /2, and the change in the in- 
terplanar distance between consecutive planes i and j 
as Adij = Pi — Pj. The rumpling parameter of layer i 
describes the movement of the ions with respect to the 
mean position of each atomic plane and corresponds to 
Vi = [$z (Mi) — S z (Oi)] /2. It is positive when the cation 
Mi is above the oxygen, and negative otherwise. 

Fig. [l] shows a schematic view of the atomic relaxations 
for both Ba and Sr-interfaces, when the in-plane lattice 



IV. ELECTRONIC STRUCTURE AT THE 
INTERFACE 

In Figure we report the energy band structure along 
a selected high symmetry line in the first Brillouin zone 
for the bulk alkaline-earth oxides (rX-line) and cubic 
bulk perovskite structures (rR-line). Only bands close 
to the Fermi level are represented. The valence bands 
are mainly composed of O 2p states that, in the case of 
the perovskites, show significant hybridizations with Ti 
3d orbitals. 

All the alkaline-earth and perovskite oxides we con- 
sider are insulators (experimental gaps above 3 eV). 
Within the cubic symmetry (in the absence of strains) 
and neglecting spin-orbit couplings, the computed band 
gap for both BaTi0 3 and SrTi0 3 is indirect (R -> T). 
The top-most valence and the bottom-most conduction 
bands are three times degenerated at these high symme- 
try points. Under the same conditions, the gap is also 
indirect in SrO, with three degenerated upper most va- 
lence bands at T and a single lowest conduction band at 
X, whereas BaO exhibits a direct gap, between single 
bands at X. 

In Table [v| we report experimental and theoretical 
bands gaps (within LDA) for all the materials involved 
in our study. We see that, due to the well known DFT 
"band gap problem" , the theoretical values are understi- 
mated by about 50 % in each case. Nevertheless, it is 
usually accepted that this error can be roughly compen- 
sated by an appropriate shift of the conduction bands 
which should not affect the conclusions of the character 
of the gap reported in this Section. 

A uniaxial strain along (001) lowers the symmetry of 
the perovskites from Pm3m to P4mm. This translates 
into a splitting of the top of the valence bands into a 
singlet and a doublet. The singlet is above (below) for 
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+0.017 
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-0.008 

-i - 
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-i- O (3) ri = -0.015 

+0.012 

4_ ! O (2) T| = +0.004 

+0.009 

0(1) T| = 



symmetry AO plane 



FIG. 1: Schematic view of the atomic relaxation for the bottom half of BaO/BaTi03 (panel a), and SrO/SrTi03 (panel b) 
supercells. Dashed lines correspond to the reference positions of the atomic planes, and the full lines are the mean position in 
the relaxed structure. Changes in the interplanar distance are written in A. The atoms (A or Ti, depending on the layer, at 
the left and O at the right) move in the direction indicated by the arrow. The rumpling parameter, r), is expressed in A. The 
size of the heterostructure corresponds to n~6, m=5. 
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FIG. 2: Bulk band structures of cubic BaO, BaTiOs, SrO, 
and SrTi03 at the theoretical lattice constant. The zero of 
energy has been asigned to the top of the valence band in 
each case. Only the bands closed to the gap are plotted. 



a compressive (tensile) strain. For the alkaline-earth ox- 
ides, the symmetry reduces from Fm3m to I4/mmm. The 
top of the valence band of SrO is therefore split but, in 
this case, the doublet is above (below) the singlet for a 
compressive (tensile) strain. Spin-orbit couplings (not 
considered in this work) might introduce further split- 
tings. 

Fig. | shows the Projected Density Of States (PDOS) 
on the different atoms (sum of the projections of the DOS 
on all the atomic orbitals of the given atom) as a func- 



TABLE V: Theoretical (Ef*°) and experimental {E^f) 
band gaps in eV for the materials involved in our simulations. 
The theoretical value, within LDA, has been calculated at the 
theoretical lattice constant. 





BaO 


BaTiOs 


SrO 


SrTiOs 


r^theo 
a gap 


1.75 


1.63 


3.03 


1.67 


rpexpt 
^gap 


4.8 


3.2 6 


5.7 " 


3.3 c 



"W. H. Strehlow, and E. L. Cook, Ref. Q 
b S. H. Wemple, Ref. El 

C R. A. McKee, F. J. Walker, and M. F. Chisholm, Ref. H 



tion of the depth of the layer inside the material for the 
BaO/BaTi0 3 and SrO/SrTi0 3 interfaces. The main con- 
clusions that can be extracted are as follows: (i) the ab- 
sence of any interface induced gap states clearly demon- 
strates the semiconductor character of the heterostruc- 
tures; (ii) the features of the PDOS on the alkaline-earth 
and the O atom at the interface layer (labelled as 4) are 
much closer to the ones displayed in bulk-AB03 than 
in bulk-AO, showing a dominant ABO3 character of the 
interface; (iii) the PDOS converges very quickly to the 
bulk properties and many of the bulk features can be re- 
covered even at the atomic layers closest to the interface; 
(iv) atomic relaxations have small effects on the shape of 
the PDOS, as can be seen comparing the solid and dot- 
ted curves in the figure. Only a shift in the SrO layers 
towards the zero energy (chosen as the top of the valence 
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FIG. 3: Projected Density Of States on all the atoms as a function of the distance to the interface for the BaO/BaTi03 (panels 
a and b) and SrO/SrTi03 (panels c and d) heterostructures. Full lines represent the projection for the relaxed geometry and 
dotted lines for unrelaxed coordinates. Atomic layers are identified as in Fig. jjj. Projected Density Of States of bulk AO and 
ABO3 are also plotted for comparison. All the energies have been rigidly displaced in order to align the top of the valence band 
(vertical line) with zero. The imposed ay was set up to the theoretical one of Si (5.389 A). 



band in each case) is noticeable. This effect is a direct 
consequence of the relaxation-induced interface dipole 
discussed in Section H. The different magnitude of the 
dipole between the Ba and Sr-interfaces explains why the 
shift is almost negligeable in the Ba-heterostructure. 



V. BAND OFFSET 



One of the most important physical quantities that 
characterize the interface between semiconductors or in- 
sulators is the band offset, i. c., the relative position 
of the energy levels on both sides of the interface. The 
valence-band offset, VBO (respectively conduction-band 
offset, CBO) is defined as the difference between the po- 
sitions of the top of the valence bands (respectively the 
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bottom of the conduction bands) of the two materials. 
These band discontinuities play a fundamental role in 
calculating the transport properties through heterojunc- 
tion devices. 

The determination of these offsets from first-principles 
cannot be achieved from a direct comparison of the corre- 
sponding band edges in the two compounds as obtained 
from two independent bulk band-structure calculations. 
The reason is the lack of an intrinsic energy scale to re- 
fer all the energies: in a first-principles simulation, the 
hamiltonian eigenvalues are referred to an average of the 
electpastatic potential that is ill-defined for infinite sys- 
temsEJ (it is only defined to within an arbitrary con- 
stant). Consequently, together with the eigenvalue dif- 
ference, we must consider the lineup of this average be- 
tween the two materials. This potential shift depends on 
the dipole induced by the electronic charge transferred 
from one part of the interface to the other after the in- 
tcrfacial hybridization (the electronic charge density of 
each system will decay into the other in an, in principle, 
unknown way). The transfer of charge depends not only 
on the materials that constitute the interface, but also 
on the particular orientation, so the lineup can only be 
obtained from a self-consistent calculation on a supercell 
including both materials. 

Therefore, from the theoretical point of view-rfche band 
offsets (BO) are usually split into two terms e3c3: 

BO = AE VjC + AV (2) 

The first contribution, AE V (resp. AE C ), is referred 
to as the band- structure term. It is defined as the dif- 
ference between the top (resp. bottom) of the valence 
(resp. conduction) bands as obtained from two inde- 
pendent standard bulk band-structure calculations at the 
same strained geometries as in the supercell construction. 
Within LDA, only a first estimate of the band-structure 
term can be obtained, AE^® A . To get more accurate 
results, a correction dealing with many-body effects in 
the quasiparticle spectra should be added: 



Knowing the relative position of the valence bands, we 
simply add the experimental band gaps (see Table |v|) 
to obtain the discontinuities for the conduction bands 
{AE c = AE^ DA + AE^ a f). 

The second term, AV, is the lineup of the average 
of the electrostatic potential through the heterojunction. 
This macroscopic quantity summarizes all the intrinsic 
interface effects, such as the chemical composition, struc- 
tural details and orientation. To obtain it, we start from 
the total (ionic plus electronic) microscopic electrostatic 
Hartree potential, output of the self-consistent supercell 
calculation (in this Section we will define the zero-enepgy 
level as the average of this potential in the unit cpll P.t) .. 
Then, we apply the double-macroscopic averagecIUj'Ej 
technique. It consists of performing first the average of 
the electrostatic potential over planes parallel to the in- 
terface, and then averaging the obtained quasiperiodic 
one-dimensional function with two step-like filter func- 
tions whose lenghts, l\ and I2, are determined by the 
periodicity of the constituents. Here l\ and I2 have been 
set up to the distance between equivalent AO and Ti02 
planes in the alkaline-earth oxide and in the perovskite 
respectively. A full description of the method to the 
AO/ABO3 heterostructures can be found in Ref. |54|. 
The resulting profile of the macroscopic potential is flat 
on both sides far enough from the interface {bulk-like re- 
gions). AV is defined as the difference between these 
two plateau values (see Fig. [i]). The lineup should be 
independent of the length used in the filter functions. 
However we have checked how doubling the size of the 
step-like functions introduce a numerical uncertainty in 
AV of the order of 30 meV. This is the main source of 
inaccuracy in our calculations of the band offsets. 

It is worth noticing that neither AE V nor AV have 
any physical meaning by their own, being pseudopoten- 
tial dependent numbers. Only the sum of both is phys- 
ically significant and-quite independent of the choice of 
the pseudopotential E3. 



A. First principles results 



AE V , C = AE V c 



LDA 



AEf 



(3) 



Self-energies corr.eGtions are often obtained within the 
GW approximationea. They strongly modify the descrip- 
tion of the conduction bands, and tend to solve the "band 



gap problem" mentioned in Section IV . Even the valence 
band energies might be subject of certain errors, specially 
in oxidesti3. Unfortunately, no accurate GW data are 
currently available for AO and ABO3 compounds. Only 
model GW calculations have been performed recently for 
SrO and SrTi03 and with limited successEj. To overcome 
the problem, we make the approximation that the errors 
in the valence bands are smaller than those for the con- 
duction bands and of the same order of magnitude for 
the two compounds taking part in the heterostructures 
so that they tend to cancel each other (AE^ orr = 0). 



Fig. H shows a schematic representation of the 
band structure discontinuities for the BaO/BaTi03 and 
SrO/SrTi03 heterostructures, both of them calculated 
fixing the in-plane lattice constant at the theoretical one 
of Si. Band splittings steming from strain are taken into 
account in the figure. We use the same sign convention 
as Van de Walle and Martin in Ref. |55|: a positive value 
of the band offset for the discontinuity at a junction A/B 
corresponds to an upward step in going from A to B. 

From the figure, we conclude that Ba-interface is 
type-II, with both the valence and conduction bands of 
BaTi03 falling in energy below the corresponding ones 
of BaO. Within the accuracy of our calculations, BaO 
and BaTi03 topmost valence bands are almost aligned 
(an offset of only -0.06 eV is predicted), so the barrier in 
the conduction bands is mainly due to the difference in 
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FIG. 4: Schematic representation of the valence-band offset (VBO) and the conduction-band offset (CBO) for BaO/BaTi03 
(panel a), and SrO/SrTi03 interface (panel b). E„, E c , and E^p stand for the top of the valence band, the bottom of the 
conduction band and the experimental band gap respectively. Values for E„, measured with respect to the average of the 
electrostatic potential in each material, are indicated. The solid curve represents the profile of the macroscopic average of the 
total electrostatic potential across the interface. AV stands for the resulting lineup. The in-plane lattice constant was set up 
to the theoretical one of Si (5.389 A). The size of the supercell corresponds to n = 6 and m = 5. 



the band gaps of both materials and, inferred from the 
experimental values, amounts to -1.62 eV. 

Sr-interface is type-I, meaning that the band gap of 
SrTiC>3 lies completely inside the gap of SrO. An upward 
step of +0.18 eV for the valence bands is theoretically 
expected, which implies a CBO of -2.22 eV. 

A rougluestimate of the valence band offset was already 
accessible E3 by identifying in Fig. [3] the position of the 
top of the valence band in the PDOS for the O atom 
at both symmetry planes, in the bulk-like regions of the 
materials that constitute the interface (numbered as 1 for 
AO and 9 for ABO3). The values deduced from the figure 
are -0.09 eV for the Ba-interface (BaO above) and +0.28 
eV for the Sr-interface (SrTiOa above), close to those ob- 
tained using the macroscopic average technique. However 
these numbers must be taken with carecS: this method to 
compute band offsets requires calculations with a higher 
number of special fc-points than those needed to converge 
the charge density or the potential lineup. 

It is important to note here the crucial role played by 
the atomic relaxations at these polar interfaces. As was 
pointed out in section [II, after the relaxation process 
an extra dipole appears at the junction that, modifies the 
electrostatic lineup accross the interface E3 and, conse- 
quently, the band offsets : 



S(AV) - 



7 *(T) 



Au Kl n 



On " e c 

KOi 



(4) 



where 8 (AV) is the change in the electrostatic lineup 

*(T) 

along z-direction due to the atomic displacements, Z K a p 



TABLE VI: Valence-band offsets (VBO) for BaO/BaTi0 3 
and SrO/SrTi03 interfaces. Values are reported at different 
in-plane lattice constants, ay. AE V and AV stand for, re- 
spectively, the band structure term and the line up of the 
electrostatic potential contributions to VBO. The size of the 
heterostructures corresponds to n = 6, m = 5. 



BaO/BaTi0 3 


SrO/SrTi0 3 


an (A) 5.389 5.430 5.583 5.665 


5.389 5.430 5.522 


AV -0.834 -0.833 -0.746 -0.755 
AE V 0.772 0.807 0.600 0.560 
VBO -0.062 -0.026 -0.146 -0.195 


-1.217 -1.190 -1.128 
1.401 1.340 1.209 
0.184 0.150 0.081 



is the Born effective tensor of atom k, Au Ka its displace- 
ment along cartesian direction a during the relaxation 
and Coo the optical dielectric constant. Looking at the 
magnitude of the atomic displacements, it is reasonable 
that the change should be more remarkable for the Sr- 
interface than for the Ba one. From our ab-initio calcu- 
lations, and for the same supercell used to get results in 
Fig. |[ we observe a change in AV of -0.67 eV for the 
Sr-interface (from 1.16 eV for the unrelaxed geometry to 
0.49 eV after the relaxation), whereas in the Ba-interface 
the deviation amounts to -0.11 eV (from 0.44 eV to 0.33 
eV). This emphasizes the importance of performing ac- 
curate first-principles atomic relaxations for correct pre- 
dictions of the barriers. 

To what extent do these discontinuities change with 
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the in-plane lattice constant? This is an important ques- 
tion because a dependence with strain would allow us to 
tune the band offsets (for example, replacing the Si sub- 
strate by Geo in order to impose a different lattice param- 
eter throughout the interface) depending on the required 
values for a given device. In order to check this point, 
we have carried out calculations at different in-plane lat- 
tice constants. In Table VI and Fig. || we summa- 
rize the results for both, Ba and Sr-heterostructures. In 
both cases a variation by about 0.1 eV in VBO with the 
in-plane lattice constant is observed, mainly due to the 
band-structure term (consequence of the strain-induced 
splittings of the top valence-band manifold), as it hap- 
pens for other lattice- mismatched, isovalent, common an- 
ion interfaces E3. The change is almost linear, and tends 
to lower the energy of the valence bands of the ABO3 
perovskite with respect to the AO alkaline-earth oxide. 

The band-structure term displays a linear behaviour 
with strain for the Sr-interface (see Fig. ||-panel a). 
The anomalous behaviour of AE V for the Ba-interface is 
due to a modification in the character of the top of the 
valence band of BaO under strain. It changes from X 
when BaO is compressed to Z when it is expanded. This 
transformation occurs for a lattice constant around 5.43 
A (theoretical lattice parameter of BaO). In Fig. ||(a) we 
plot the difference between the top of the valence band 
of BaTi03 and the highest occupied state at X and Z 
of BaO. The crossing point is clearly identified in the 
figure. No extra changes in the linear beaviour of AE V 
are expected for longer lattice constants. 

The almost-linear change in the lineup term can be ex- 
plained according to an analytic scaling law proposed in 
Rcf. |5^. Once AV is known for a reference configura- 
tion with an in-plane lattice constant ay, then, supposing 
an uniform strain throughout the structure, AV for any 
other strained configuration a H can be extrapolated from: 



AV 



1 



AV 



(5) 



where e\\ = (a,n/a\\ — 1). Fig. |^(b) shows a comparison 
of the first-principles and extrapolated values, where the 
heterostructure at the in-plane lattice constant of Si has 
been chosen as the reference configuration. Results are 
in good agreement within the numerical accuracy of the 
ab-initio results. 

In summary, the VBO varies almost linearly for a large 
range of in-plane strains. The only deviation is observed 
for BaO/BaTi03 and is explained by a change of char- 
acter of the BaO gap under compression. 



VI. INTERFACE WITH SI 

As it was pointed out in the Introduction, AO/ABO3 
interface is only a part of the gate stack of technologi- 
cal interest for the semiconductor industry. AO acts as 
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FIG. 5: Dependence with in-plane strain of the valence- 
band offset (VBO), and of its decomposition into the band- 
structure (A_B„), and lineup (AV) terms. Squares and cir- 
cles represent, respectively, the first-principles results for the 
SrO/SrTi03 and BaO/BaTiOa interfaces. In panel (a), grey- 
filled (respectively black-filled) circles stand for the diffrence 
between the top of the valence band of BaTiOs and the high- 
est occupied state at X (respectively Z) point in BaO. Lines in 
panels (a) and (c) (dashed for Sr and full for the Ba-interface) , 
are a guide to the eye. Lines in panel (b) represent the results 
of the anayltic scaling law proposed in Ref. E9I 



a buffer layer between the Si substrate and the high- 
k perovskite. The whole heterostructure, epitaxially 
grown following the McKee- Walker process H is made of 
Si/ASi 2 /AO/AB0 3 As it will be emphasized in this 
Section, the role of the buffer layer is not only the passi- 
vation of the Si substrate, but also the efficient tuning of 
the offsets between the perovskite and the channel. 

Combining our results with various data available in 
the literature, we can estimate the band discontinuities 
along the whole heterostructures of technological inter- 
est as summarized in Fig. ||. Previous theoretical cal- 
culations of the band offsets between the alkaline-earth 
oxide AQ and Si have been reported recently (Si/BaOEil, 
Si/SrOE2l. In addition, we can find in the literature theo- 
retical estimations for the Schottky barriers between per- 
ovskites and prototypical metallic electrodes (SrTiOs/Pt 
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0, or BaTiOa/SrRuojl). 

We must notice that, altough most of the previous 
works (except Ref. [n]) have been done using the same 
basic approaches (DFT, pseudopotentials, fc-point sam- 
plings, supercells, etc), they differ in the details of the 
calculations (for example, they have been carried out at 
different in-plane lattice constant, and the size of the 
supercells, or the parameters used to generate the pseu- 
dopotentials might change from one to the other). So, 
only a rough estimate of the barriers can be deduced 
from the comparison and any quantitative conclusion is 
beyond the scope of this Section. 

Within LDA_plus GW corrections, Boungiorno- 
Nardelli et al. E3 investigated the structural and elec- 
tronic properties of the Si/SrSi2/SrO interface. They 
predicted a VBO between Si and SrO of -1.91 eV for the 
most stable interface configuration. Using the experi- 
mental gaps to locate the conduction bands, it translates 
in a CBO of 2.69 eV. The Schottky barrier <\> n (difference 
between the Fermi level and the bottom of the conduc- 
n band) between SrTiOa and Pt has been evaluated 
to -0.89 eV, which implies a barrier height 4> p (differ- 
ence between the Fermi level and the top of the valence 
band) of 2.41 eV. These results are summarized in the 
first panel of Figure ||. 

Through first-principles_gradient-density-functional 
calculations Gulleri et al. l£J focused on the structure, 
energetic and band offsets of the Si/BaO interface. For 
the favoured geometry, they obtained a VBO of -2.0 
cV. Some of us evaluated the Schottky barriers between 
BaTiOa and SrRu03 (a typical metallic electrode in fer- 
roelectric devices E5) to be equal to <j) p = +1.44 cV and <ft n 
— -1.76 eV. These results are summarized in the second 
panel of Figure 0. 



For both stacks, we can clearly see how the problem of 
the large carrier injection (expected for the perovskite in 
direct contact with SitiJ) is overcome by the use of the AO 
layer. The electrostatic barriers for both the electrons 
and holes, between the gate electrode and the channel 
are large enough to prevent carrier injections and to push 
the use of ABO3 perovskites to a prominent position to 
replace silica as the gate dielectric oxide in MOSFETs. 

VII. CONCLUSIONS 

We have studied structural and electronic properties 
of BaO/BaTi0 3 and SrO/SrTi0 3 interfaces from first- 
principles. Atomic relaxations have been performed. In- 
terface dipoles, due to the opposite motion of the anion 
and cation atoms at the interface, appear for both het- 
erostructures. No interface electronic states are induced 
in the band gap. The character of the AO layer at the in- 
terface is mainly perovskite-like. Under the experimental 
strain conditions, the valence bands of BaO and BaTiOs 
are almost aligned (within the accuracy of our calcula- 
tions), whereas a slightly larger barrier is predicted for 
SrO/SrTi03. Interface dipoles, induced by atomic relax- 
ations, have a strong effect on the band alignments at 
the interface. A nearly linear variation of the VBO with 
in-plane strain is observed. 

Gathering together our results and various data 
available in the literature, we make a guess for 
the band alignment of whole Si/SrO/SrTiOs/Pt and 
Si/BaO/BaTiOs/SrRuOa structures. In both cases large 
enough electrostatic barriers for electrons and holes be- 
tween the gate electrode and the channel are estimated, 
preventing the injection of carriers and suggesting that 
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both perovskites compounds are promising candidates to 
replace silica in MOSFETs. Our results should be con- 
firmed by more accurate calculations for the whole het- 
crostructure. 
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